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^ , Abstract 

, The equations for spin evolution of a particle in a storage ring are obtained considering con- 

' tributions from the tensor electric and magnetic polarizabilities of the particle along with the 

contributions from spin rotation and birefringence effect in polarized matter of an internal tar- 
T-H ' get. Study of the spin rotation and birefringence effects for a particle in a high energy storage 

, ring provides for measurement both the spin-dependent real part of the coherent elastic zero-angle 

' scattering amplitude and tensor electric (magnetic) polarizabilities. 

O ■ 1 INTRODUCTION 
\0 . 

I Investigation of spin-dependent interactions of elementary particles at high energies is a very important 
• part of programs for scientific research at storage rings T', '2'. Such studies are being carried with the 
use of polarized beams and polarized targets. Dependence of scattering cross-sections on the particle 
OhI spin is the subject of much studies. The experiments for measuring the spin-dependent part of the 
forward scattering amplitude are in preparation now fT.. 

It should be mentioned that it is well known in experimental particle physics how to measure a 
. ^ , spin-dependent cross-section. However, measuring of the spin-dependent part of the forward scattering 
^ I amplitude is the complicated challenge. 

' It was shown in [Sl-|n] that there is an unambiguous method, which makes the direct measurement 

of the real part of the spin-dependent forward scattering amplitude in the high energy range possible. 
This technique is based on the effect proton (deuteron, antiproton) beam spin rotation in a polarized 
nuclear target and on the deuteron birefringence effect i.e. phenomenon of deuteron spin rotation 
and oscillation in a nonpolarized target. This technique is based on measurement of the angle of spin 
rotation of a high energy proton (deuteron, antiproton) in conditions of transmission experiment. 

The similar phenomenon for thermal neutrons was theoretically predicted in |li]j and experimen- 
tally observed in jllj-|13j (the phenomena of nuclear precession of neutron spin in a nuclear pseudo- 
magnetic field of a polarized target). 

When considering particles moving in a storage ring one should take into account infiuence of 
electromagnetic fields, which exist in the storage ring, on behavior of the particle spin. 

According to |14[ I15j . particle spin behavior in electromagnetic fields can be described by the 
Bargmann-Michel-Telegdi equation. When the particle has the non-zero electric dipole moment 
Bargmann-Michel-Telegdi equation should be supplied with additional terms, which describe inter- 
action of the EDM with the electromagnetic field in the storage ring. This interaction is much weaker 
than conventional interaction of the particle spin with the electromagnetic field. Nevertheless, the 
modern experimental technique provides to measure such weak interaction of EDM with the elec- 
tromagnetic field |16[ I17j and to carry out experimental search of the deuteron EDM at the level 
d ~ 10~^^e • cm. 
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In the present paper it is shown that modern experimental technique ahows to study particle spin 
rotation in a polarized internal target of a storage ring and, thus, to measure the spin-dependent part 
of the forward elastic coherent scattering amplitude. It is shown that dynamics of the particle spin in 
a storage ring in the experiment with an internal target can not be completely described by the BMT 
equation. Moreover, it is shown that dynamics of the particle spin in a storage ring in the experiment 
with the accuracy necessary for the EDM search for particles with the spin S > 1 (deuteron) can not 
be completely described by the BMT equation, too. This is due to the particle possesses the electric 
and magnetic tensor polarizabilities. The equations that can be applied in this case for description of 
spin evolution are obtained and additional contributions in these equations are analyzed. 



2 Particle spin rotation in an electromagnetic field in a storage ring 



Thus, let us consider now a particle with the spin S moving in the electromagnetic field of a storage 
ring. Interaction of the particle magnetic moment with the electromagnetic field presenting in the 
storage ring leads to rotation of the particle spin with respect to the momentum direction. To describe 
particle spin evolution it is used the Bargmann-Michel-Telegdi equation |14[ II 5j as follows: 
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m is the mass of the particle, e is its charge, P is the spin polarization vector in the deuteron rest 
frame, 7 is the Lorentz- factor, P = v/c, v is the particle velocity, a = {g ~ 2)/2, g is the gyromagnetic 
ratio, E and B are the electric and magnetic fields in the point of particle location. 

The parameter a for protons is o = 1.79 and for deuterons a = —0.14. As a result, for the typical 
field in a storage ring {B « 10^ ^ 10^ gauss) the spin rotation frequency is O ~ 10^ -i- 10^ sec ^. 

If a particle possesses an intrinsic dipole moment then the additional term that describes spin 
rotation induced by the EDM should be added to ^Hl 
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where d is the electric dipole moment of a particle. 

As a result, evolution of the deuteron spin due to the magnetic and electric dipole momenta can 
be described by the following equation: 



dP_ 
'dt 



e 

mc 



P X 



a + -]B 

7. 



7 



7 + 1 



P-B)P 



7 



7 + 1 



PxE 



+ d 



Px (c(3x B + E) . (4) 



The typical frequency uJd of precession caused by particle EDM interaction with the electromagnetic 
field is many orders less than 0. In conditions of the experiments |16[ I17j. which are planned for the 
deuteron EDM search at the level d ~ lO^^^e • cm, the precession frequency is ~ 10~^ • 10~^ sec~^. 
Nevertheless, according to |16| I17j. the change of deuteron spin direction in the storage ring can be 
measured even when caused by such small uJd- 

When considering particles with the spin S > 1, the more so since the accuracy of the experiments 
is expected to be high, the particular attention should be focused on behavior of the particle spin in 
the storage ring, because it is strongly influenced by several additional interactions [HlElElEni- Note 
that a particle with the spin S > 1 has the electric and magnetic polarizabilities. Particle interaction 
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with the electromagnetic field in the storage ring aroused by its tensor electric polarizability can 
be expressed as follows (suppose that E and B are orthogonal to the particle momentum - this case 
is realized in the storage ring): 

= -2'^ik{Eeff)i{Eeff)k, (5) 

where Ojfc is the electric polarizability tensor of the particle, -Ee// = {E + (3 x B) is the effective electric 
field in the point of particle location. The expression can be rewritten as follows: 

= asElff - arEljf (Sue) , ue = (6) 
^ ^ \E + P X B\ 

where as is the scalar electric polarizability and is the tensor electric polarizability of the particle. 

The mentioned interaction results in rotation and oscillations of the spin of a particle moving in 
an electric field [HI El El HOI- The frequency of this rotation ujte is determined by the tensor electric 
polarizability ot- Theoretical evaluations \21\ give for the deuteron polarizability ax ~ 10"^'' cm^. 
Therefore, for a deuteron in a storage ring the typical energy (in the frequency units) of spin interaction 
with the electric field due to the tensor electric polarizability is about ojte ~ 10"^ in the field 

A particle with the spin 5 > 1 also has the magnetic polarizability, which is described by the 
magnetic polarizability tensor and interaction of the particle with the magnetic field due to the 
tensor magnetic polarizability is as follows: 

= --^PikiBeff)i{Beff)k, (7) 

where {Beff)i are the components of the effective magnetic field ^e// = {B — (3 x E); the energy of 
interaction Vg (|7|) can be expressed as: 

= PsBlf^ - hBlf^ (Sub) \nB= (8) 
V / \B - j3x E\ 

where [is is the scalar magnetic polarizability and /3t is the tensor magnetic polarizability of the 
particle. 

Similarly interaction Vg, interaction Vg arouses rotation and oscillations of the particle spin with 
the frequency Uj,. According to calculations 21 the tensor magnetic polarizability can be evaluated as 
/3t ~ 2-10"^'^ cm?. Therefore, this contribution provides rotation with the frequency ~ 10~^ sec~^ 
in the field -Be// ~ 10^ gauss. 

As it is obvious the frequencies lote and ujj, are much higher than the frequency u;^, caused by 
the deuteron EDM (if it is nonzero). Therefore, the influence of the tensor polarizabilities on particle 
spin evolution in a storage ring should be considered with particular attention. 



3 The index of refraction and effective potential energy of particles 
in medium. 

Recall now that a storage ring can contain matter inside (gas, jet target). Presence of such a target 
influences on spin behavior j^l 1181 1191 l2Uj along with the influence from the electromagnetic fields. 
To make understanding of target influence on spin behavior clear let us recollect the following. Close 
connection between the coherent elastic scattering amplitude at zero angle /(O) and the refraction 
index of medium n has been established as a result of numerous studies (see, for example, |22| I23j): 

27rA^ . . ^ 

n = l + -r^/(0) (9) 
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where is the number of particles per cm^, k is the wave number of the particle incident on the 
target. 

The expression (0 was derived in assumption that n — 1^1. If/c— >0 then (n — 1) grows and 
expression for n has the form 

n2 = l + ^/(0) 

Let us consider particle refraction on the vacuum-medium boundary. 

The wave number of the particle in vacuum is denoted k. The wave number of the particle in 
medium is k' = kn. As it is evident the particle momentum in vacuum p = hk is not equal to the 
particle momentum in medium. Therefore, the particle energy in vacuum + m^c^ is not 

equal to the particle energy in medium Emed = y/h^k'^n^c^ + m^d^ 

boundary vacuum-medium 



vacuum 

E 



medium 
^med 



Figure 1: The particle energy in vacuum E is not equal to that in medium Emed- 

The energy conservation law immediately necessitates to suppose that the particle in medium 
possesses an effective potential energy V^ff (see the detailed theory in This energy can be easily 

found from the evident equality 

E = Emed + Veff 

i.e. 

Veff = E- Emed = Nf{E, 0) = {2TrfNT{E), (10) 

m7 

where it is used f{E,0) = -{27rf ^ T{E) = -{2^Y ^ r(^), r(^) is the T-matrix 

Above we considered the rest target. But in storage rings moving bunches can be used as a target. 
Therefore we should generalize the expressions ()9I1U|) for this case. Thus, let us consider the collision 
of two bunches of particles. Suppose that in the rest frame of the storage ring the particles of the first 
beam have the energy Ei and Lorentz- factor 71, whereas particles of the second beam are characterized 
by the energy E2 and Lorentz- factor 72. Let us recollect that the phase of a wave in a medium is 
Lorentz- invariant. Therefore, we can find it by the following way. Let us choose the reference frame, 
where the second beam rests. As in this frame particles of the second beam rest, then the refraction 
index can be expressed in the conventional form ©: 

n; = i + ^/(£;;,o), (11) 

where N2 = ^2^N2 is the density of the bunch 2 in its rest frame and is the density of the second 
bunch in the storage ring frame, fc^, E'^^ are the wavenumber and energy of particles of the first bunch 
in the rest frame of the bunch 2, respectively. Suppose the length of the bunch 2 in its rest frame is 
-L, then L = 72 /, where I is the length of this bunch in the storage ring frame. 

Now the change of the phase of the wave caused by the interaction of the particle 1 with the 
particles of bunch 2 can be found: 

ct> = k[{n[ - l)L = ^^f{E[,0) L = ^^f(E[,0)k[ I, (12) 
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n, = l + ^-^^m,0) (15) 



It is known E3]that the ratio "^^^J'^^ is invariant, therefore, ^^^}'^^ = ^^^['^^ i where /(i?i,0) is the 
amphtude of elastic coherent forward scattering of the particle 1 by the moving particle 2 in the rest 
frame of the storage ring. 

As a result 

^ = ^f{El,0) I = ^f{E,,0) Vrel t, (13) 

fci ki 

where Vrei is the velocity of relative motion of the particle 1 and bunch 2 (for opposing motion 
i^rel = (^^1 + ^2)(1 + ^^jS^)~^)j t is the time of interaction of the particle 1 with the bunch 2 in the rest 
frame of the storage ring. 

The particle with the velocity vi = passes the distance z = vi t over the time t. It should 

be noted that the path length z differs from the length of the bunch 2, because it moves. Expression 
(|13j) can be rewritten as: 

</)=^/(ii;i,0) ^z = fci(ni-l)z, (14) 
ki vi 

where the index of refraction of the particle 1 by the beam of moving particles 2 is: 

ki^ vi 

When V2 = ^1 the conventional expression Q follows from (|15() . 

Now the effective potential energy Veff being acquired by a particle 1 when colliding with the 
particles of bunch 2 can be found: 

Veff = Ei-Ei w = Ei- Jv\(?n\ + m\c^ = -27:h^ N2VreJ-^^^ = -27Th^ N2Vre /^^y^\ (16) 

\ 1- 1- ^ Pi p'^ 

therefore 

Veff = -^^^f{E[,0) = {2nfN2T{E[), (17) 
mi7i72 

where E[ = mic^7i72 is the energy of the particle 1 in the rest frame of the bunch 2, pi is the particle 
1 momentum in the storage ring frame, while p[ is the particle 1 momentum in the rest frame o the 
bunch 2, p'^ = it is taken into account that v'l = Vrei- 
When obtaining (|T6|) it was used |ni — 1| <C 1. 

Let us consider now what happens when the particle possesses a spin. In this case the amplitude 
of the zero-angle scattering depends on the particle spin and, as a consequence, the index of refraction 
depends on the particle spin and can be written as: 

ni = l + — ^ f{Ei,0), (18) 

ki vi 

where / (-Ei,0) = TrpjF (0), pj is the spin density matrix of the scatterers, -F (0) is the operator of 
the forward scattering amplitude, acting in the combined spin space of the particle and scatterer spin 
J. 

According to the above (see (|10I17|) ') a particle possesses some effective potential energy V in 
matter. If the amplitude /(O) of particle scattering depends on the particle spin, then the effective 
energy depends on the spin orientation 0-11]: 

"117172 

Thus, for example, the amplitude /(O) of the spin 5 = 1 particle (for example, deuteron) scattering 
even in an unpolarized target depends on the particle spin and can be written as: 

f{E[,Q) = d + di{Snf. (20) 
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where S is the deuteron spin operator, n is the unit vector along the deuteron momentum k. 

Substituting the expression for /(O) (PUj) in ((T^ one can obtain for a particle with the spin S = 1: 

Veff = N2id + dJSn) ) . (21) 

mi7i72 V ^ W 

Let the quantization axis z is directed along n and M denotes the magnetic quantum number. 
Then, for a particle in a state that is the eigenstate of the operator Sz of spin projection onto the 
z-axis, the effective potential energy can be written as: 

^eff = ^^iV2 (d + diM^) . (22) 

mi7i72 

According to (|22|) splitting of deuteron energy levels in matter is similar to splitting of atom energy 
levels in the electric field aroused from the quadratic Stark effect. Therefore, the above effect could 
be considered as caused by splitting of the spin levels of the particle in the pseudoelectric nuclear field 
of matter. 

Comparison of (|21() with provides to conclude that interactions (|21() and © are similar, there- 
fore we can observe the effect of spin rotation and oscillations for a particle with S > 1 passing through 
nonpolarized matter (birefringence effect) [HI EI- Henceforth, for the expression (|21|) we will use the 
notation Vg"""^ = Veff. 

4 Spin rotation of proton (deuteron, antiproton) in a storage ring 
with a polarized target 

Let us consider now the experiments deals with the use of polarized beams and targets. In this case 
the amplitude of the elastic coherent scattering at the zero angle depends on the vector polarization 
Pt of the target nuclei (if the spin of the target nuclei J > 1, then the addition depending on the 
target tensor polarization also appears, but this addition will be omitted in below considerations, the 
general case is considered in CI ) ■ 

For the sake of concreteness suppose the target to be rest (72 = 1, 71 = 7). 

The contribution to the amplitude /(O) proportional to the vector polarization of the target nuclei 
can be expressed as: 

/(Pt,0) = AiSPt + A2{Sn){nPt). (23) 

Contributions from weak P,T-odd interactions (see 0, |18)-|2n)) are omitted here. 

Correspondingly, the contribution to the effective potential energy of particle interaction with 
matter caused by polarization of target nucleus spins looks like: 

VeffiPt) = -^N{AiSPt + A2{Sn){nPt)). (24) 
7717 

The expression Veff{Pt) can be rewritten as: 

V^^^' = VeffiPt) = -flG = -^SG (25) 
where ^ is the particle magnetic moment, 

G = '^^^N{AiPt + A2n{nPt)). (26) 
m7/x 



6 



Recall now that the energy of interaction of the magnetic moment fl with a magnetic field B is as 
follows: 

Vmag = -pB = -^SB. (27) 

Comparing H25|) and ()27() one can easily find they are identical. Therefore, G can be interpreted 
as the effective pseudomagnetic field acting on the spin of the particle in matter with polarized nuclei. 
Pseudomagnetic field G is caused by nuclear interaction of particles with scatterers. Similar particle 
spin precession in conventional magnetic field, spin also precesses in the field G (this phenomenon is 
called nuclear precession of the particle spin). This phenomenon was described for the first time for 
slow neutrons in ^H] and observed in pT]-P^. 

5 The equations describing spin evolution of a particle in a storage 
ring 

5.1 Interactions contributing to the spin motion of a particle in a storage ring 

Thus, according to the above analysis, when considering evolution of the spin of a particle in a storage 
ring one should take into account several interactions: 

1. interactions of the magnetic and electric dipole moments with electromagnetic fields; 

2. interaction of the particle with electric fields due to the tensor electric polarizability; 

3. interaction of the particle with magnetic fields due to the tensor magnetic polarizability; 

4. interaction of the particle with the pseudoelectric nuclear field of matter; 

5. interaction of the particle with the pseudomagnetic nuclear field of polarized matter. 

The equation for the particle spin wavefunction considering all these interactions is as follows: 

ih^^ = (i?o + Vedm + Vj^ + Vg + yr^' + Vg"""') ^(0 (28) 

where ^{t) is the particle spin wavefunction, 

Hq is the Hamiltonian describing spin behavior caused by interaction of the magnetic moment 
with the electromagnetic field (equation (|28|) with the only Hq summand converts to the Bargmann- 
Michel-Telegdi equation) ; 

Vedm describes interaction of the particle EDM d with the electric field; 

Vg describes interaction of the particle with the electric field due to the tensor electric polarizability; 
is the energy of interaction of the particle with the magnetic field due to the tensor magnetic 
polarizability; 

ynuci (igsCTiijgg the effective potential energy of particle interaction with the pseudoelectric field 
of the target ([^1) . 

ynucl (^ggcril^gg the effective potential energy of particle interaction with the pseudomagnetic field 
of the target (^5]) . 

5.2 Deuteron spin rotation in electromagnetic fields of a storage ring 

Let us consider a deuteron moving in a storage ring with low pressure of the residual gas (10^^'^ 
Torr) and without targets inside the storage ring. In this case we can omit the effects caused by the 
interactions V^"'^' and Vq^'^K According to the above analysis spin behavior of a deuteron, nevertheless, 
can not be described by the Bargmann-Michel-Telegdi equation. According to the above interactions 
Ve and Vb should be considered when describing spin evolution. The corresponding equations can be 
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written as follows [TH]-pn]: 



where 



dP 
dt 



P X n{d) 



2 I^tBIjj 



[uE X P]inE,k + nE,i[nE x P]k 
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= 17 + Od, 
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' E + f3x B 



n -- 

^d-Ti 

Pik = Sp p Qik, 

Qik — 2 ['S'j'S'fc + SkSi — gOjfcJ 



p is the spin density matrix of the particle, m is the particle mass, e is its charge, P is the spin 
polarization vector, Pj^ is the spin polarization tensor, P^x + Pyy + -Pzz = 0, 7 is the Lorentz- factor, 
P = v/c, V is the particle velocity, a = (g — 2)/2, g is the gyromagnetic ratio, E and S are the electric 
and magnetic fields in the point of particle location, E^jf = {E + f3 x B), B(,ff = [B — j3 x E), 



riE 



E+l3xB 
\E+0xB\ 



riB 



B-f3xE 
\B-0xE\ 



' ''^EA 



Pikn'E,k, n'^i = PiinBh ^r{d) are the components of the vector 



(r = 1,2,3 correspond to x,y,z, respectively). 
The equations for particle spin motion ()29() can be rewritten as follows: 



dP 
It 



[P X n{d)] + QrinE x n'^] + Q^[nB x n's], 



dP 



ik 



dt 



-{ejkrPij^r{d) + €jirPkj^r{d)) + ^T{[nE X P]inEk + nEi[nE X P]k) + 

+Q'^{[nB X PjiUBk + riBiinB x P]k) 



(30) 



where 



o-tE' 



, axE 



Suppose the external electric field in the storage ring = and the particle moves along the circle 
orbit. 

Let us now consider the equation H3U|) in the coordinate system that rotates with the frequency 
of particle velocity rotation. In such a system spin rotates with respect to the momentum with the 
frequency determined by {g — 2). The coordinate system and vectors v,E,B are shown in figure and 
denote the axes by x, y, z (or 1, 2, 3, respectively). 
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As a result, the system ()29() can be written as: 

nP2 - nt^P23, 



dPi 
dt 



where u^d = — 



^ = -OPi + (Q^ - VLt)Pi3 + uJdPs, (31) 
^ = nTPi2 - u;dP2, 

^ = -2nPi2, (32) 
= -0 (Pn - P22) - J^T^3 + uJdPis, 

^ = + ^tP2 - ^'^P2 - 0JdPl2 , (33) 

^ = -f)Pi3 + l^^^Pi - a;rf(P22 - Pss). 
Remember that Pn + P22 + P33 = and Pj^ = P^i- 

5.3 Contribution from the EDM and tensor polarizabilities to deuteron spin os- 
cillation 

Let us consider the system 1)311133^ more attentively. 

Suppose that deuteron has neither EDM no tensor polarizabilities: in this case the system (|31I33|) 
can be expressed: 

^ = f^^2, ^ = -f^A, ^ = 0, (34) 
^ = 217P12, = -2OP12, ^ = 0, (35) 

^ = -a (Pn - P22) , ^ = ^P2Z: ^ = -ms. (36) 

This is the conventional system of BMT equations that describes particle spin rotation with the 
frequency equal to {g — 2) precession frequency Q = ^B. The component P3 of vector polarization in 
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this conditions is equal to a constant = 0) along with the component P33 of tensor polarization 
(^ = 0). 

Suppose the deuteron EDM differs from zero. Then the above system of equations converts to: 

= ^ = -nPi + oJdPs, ^ = -uJdP2, (37) 

= 2nPu + 2uJdP23, ^ = -2nPu, ^ = -2u;dP23, (38) 

^ = -n {Pn - P22) + uJdPi3, ^ = - uJdPu, ^ = -nPis - MP22 - ^33). (39) 

From (|37ll39j) it follows that presence of the nonzero EDM makes the vertical component P3 of vector 
polarization oscillating with the frequency of {g — 2) precession 0. 

According to the idea JZj these oscillations can be eliminated by modulation of the deuteron 
velocity with the frequency 0: 

V = vo + Sv sin (nt + ipf) (40) 

here v^/ is the phase of forced velocity oscillations. 

As Eeff depends on (3 = v/c, it also appears modulated: 

Eeff = E^^fj + 6Eeffsin{nt + ^f) (41) 



Therefore uJd = — is also modulated with the same frequency. This makes the product 



I = 

<^dP2 ~ sin^(r2i + Iff). Therefore, averaging this value over the period of {g — 2) precession gives 
the result time-independent (i.e. = const) and P^it) = Ps^O) + const ■ t. For better measurement 
conditions it is important to make -P3(0) = 0. This is the reason to chose particle spin laying in the 
horizontal plane. 

All the above reasoning makes ~ const, too. Therefore, P33 also linearly grows with time 
-f33(^) = -^33(0) + const ■ t. However, when the spin lays in the horizontal plane ^33(0) 7^ 0. 

It is important to note (see below the section 5.5) that choosing spin orientation corresponding to 

cos^-!? = I (cosi? = y'^) sin-d = y'^) one gets the component ^33(0) = 0, while -P3(0) / 0. 

Therefore taking making cos"!? = we can use the component P33 for EDM measurements, 
too. 

Let us consider now the contribution from the electric and magnetic tensor polarizabilities. Then 
instead of the system (|37I39|) we should consider the system (|31I33() 

Some interesting implications follow from (|31II33|) . As it was already mentioned above in the exper- 
iments for EDM search it is planned to measure growth of the vertical component of the polarization 
vector P3. 

According to (jST|) time dependence of the vertical component of the vector polarization P3 is 
described by the equation 

^ = flrPu - 0JdP2 (42) 
at 

As it can be seen from (|31() dependence of P3 on time is determined by both the EDM and tensor 
polarizability of deuteron. It is interesting that the derivative of the tensor polarization component 
P33 does not contain contributions from the tensor electric polarizability and is proportional to the 
EDM only: 

^ = -2^rfP23 (43) 

Therefore, it is important to measure the component P33, too. According to the above spin orientation 

1 

3- 



for this case is determined by the condition cos^ "d — - 
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5.4 Contribution from the tensor magnetic polarizability to deuteron spin oscil- 
lation 

Contributions to spin rotation and oscillations from the EDM and polarizabilities are small. Therefore, 
they, being analyzed, could be considered as perturbations to the full system H31II33() and the role of 
each could be studied separately. 

The system of equations taking into account the contribution from the tensor magnetic polariz- 
ability Pt is as follows: 

!i^ = nP2- n>^P23, ^ = -nPi + nt^Pn, ^ = nP23 - n'^P2, ^ = -nPi^ + n'^Pi m 

Introducing new variables P+ = Pi + iP2 and = P13 + iP23 and recomposing equations (|1H) to 
determine P+ and G+ one obtains: 

^ = -inp+ + ini^G+, 



dt 

or 



^ = -inG+ + in'^p+, 



Let us search P+, G+ ~ e*'^*, then (|^^ transforms as follows: 



p_^ — r^p_i_ — 



uG+ = VtG+ - n'^p^ 



The solution of this system can be easily found: 

{uj - - ni^n'^ = (45) 

that finally gives 

wi,2 = n± ^Jq^^. (46) 

The solution can be rewritten as: 

P+(t) = cie-*'^i* + C2e-''^2* = |ci|e-^("^i*-'^i) + |c2|e-^(^2i-<S2) (47) 

Therefore, 

Pi{t) = \ci\ cos{LOit — 5l) + \C2\ COs{Lx)2t — 82) (48) 

P2{t) = -|ci| sin(a;it - 5i) - \c2\ sin(u;2i - ^2) (49) 

According to (|48I49() the nonzero deuteron tensor magnetic polarizability makes the spin rotating 
with two frequencies toi and W2 instead of and, therefore, experiences beating with the frequency 

/\u = uJi-u^2 = 2^0^ = 

Let us recall now that EDM interaction with the electric field makes the deuteron spin rotating 
around the direction of this field and leads to appearance of P3 component proportional to P2{t)- 

^ ~ -c^dP2. (50) 
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Therefore, 



= iOd (|ci| sin(a;it - 6i) + |c2| sin(w2i - S2)) ■ (51) 

According to the idea J3 measure the EDM one should modulate the particle velocity 
{v = vq + 5v sin {^Ift + iff)) with the frequency close to the frequency Cl oi {g — 2) precession. 

If the magnetic polarizability is equal to zero, then uji = uj2 = ^ and spin rotates in the horizontal 
plane with the frequency fi. In this case velocity modulation with the same frequency Qf = Q gives: 

dP 

^ ~ sm\nt) (52) 

and the vertical component P3 = ^^dt' linearly grows with time. 

However, coi 7^ UJ2 and velocity modulation, for example, with the frequency Q = ui provides for 
slow spin oscillation with the frequency uji — 0J2 instead of linear growth. 

According to evaluation |2J the tensor magnetic polarizability I3t 'i- 10"^*^, therefore in the field 
B ~ 10^ gauss the beating frequency Ai^ ~ 10~^. 

Measurement of the frequency of this beating makes possible to measure the tensor magnetic 
polarizability of the deuteron (nuclei). 

Thus, due to presence of the tensor magnetic polarizability the horizontal component of spin rotates 
around B with two frequencies wi, uj2 instead of expected rotation with the frequency Vl. 

This is the reason for the component P3 to experience the similar oscillations caused by the EDM. 
Therefore, particle velocity modulation with the frequency provides for eliminating of oscillation 
with the frequency Jl, but P3 oscillations with the frequency Aw rest (similarly P33). Study of these 
oscillations is necessary because they can distort the EDM measurements. 

5.5 Contribution from the tensor electric polarizability to deuteron spin oscilla- 
tion 



Let us consider now contribution caused by the tensor electric polarizability. From the system (|32|) it 
follows 



(53) 

Thus we have the equation 

^+-iVi2 = (54) 



where 0J12 = 1/4^2 -)_ Q^o;^ 2f2, because VLx^'rp <C . 
The solution for this equation can be found in the form: 

P12 = ci cos ijJi2t + C2 sin bJi2t (55) 

Let us find coefficients ci and C2: when t = the equation gives ci = Pi2(0). The coefficient 
C2 can be found from 

^^(t^0)=Wl2C2, (56) 

dt 

therefore 

C2 = T^(*^0). (57) 

UJ12 dt 
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From the equation 

dPi2 
dt 

it follows that 



{t^0) = -n {Pu{t ^ 0) - P22{t ^ 0)) , (58) 



C2 = ^ , (59) 

and 

-^11 — -^22 

Pl2 = Pl2 (0) COS L<Jl2t ^ sin 1^12^ (60) 

As a result, one can write the following equation for the vertical component of the spin P3: 

^ = nrPuit) = nT[Pi2{0)cos2nt - ^^l^^l^l^^^ sm2Qt] (61) 

As it can be seen the vertical component of the spin oscillates with the frequency 20. 
The equation describing contribution from the tensor electric polarizability and EDM to P3 looks 
like dSII): 

^ = nrPu - 0JdP2. (62) 

As P2 oscillates with the frequency 0, then the product UdP2 contains the non-oscillating terms (see 
section 5.3) and contribution to P3, which is caused by the EDM, linearly grows with time, when 
flf = 0,. If ^ 0,, then contribution to P3 caused by EDM slowly oscillates with the frequency 

Uf - n. 

It is important that modulation of the velocity v = vo + 5v sin {^Ijt + ipj) results in E(,f f oscillation 
(see (|41() ) and, therefore, E^ff also oscillates with time and appears proportional to sin^ {Qft + iff). 
As a result, the contribution to P3 caused by the tensor electric polarizability can be expressed as: 

^h^ll ^ AQt sin^ {nft + 99/) [P12 (0) cos 2nt - ^^i^W "^22(0) ^.^ ^Ot] , (63) 
where MIt = -\^{bE^f^f i.e. 

^^pi ~ -^A^TCOs (2%t + 2(^/)[Pi2(0)cos2Ot - ~ '^^^^"^ sin2r?t] (64) 

According to ()64() for a partially polarized deuteron beam the derivative depends on the 
deuteron polarization components P12 and i^iMz^. 

For simplicity let us consider a deuteron beam in the pure polarization state. In this case the 
components P12 and -^^i^o)"-^^^ ^an be written using the explicit expression for the spin wavef unctions. 
Suppose n('!?, (p) is the unit vector directed along the deuteron spin {d and are the polar and azimuth 
angles (see FigEJ). So the spin wavefunction that describes the deuteron spin state with the magnetic 
quantum number m = 1 can be expressed as follows (in the Cartesian basis): 

ax \ ^ I cos cos if — i sin ip 

Xi ("!?,¥')= I ay = — cos sin (/? + i cos 99 | (65) 

a^ ) ^2 \^ _g-j^^ 

Polarization vector can be written as: 

P = {S)= xtSxi = i[d X a*] (66) 

and components of polarization tensor 

3 2 

{Pik) = xtQikXi = -^{^ial + aka* - -}, (67) 
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where Qik is the spin tensor of rank two {Qik = ^{SiSk + SkSi — ^6ikS)). Therefore, 

3 

Pi2 = -sm2ipsm^^, (68) 
-COS Zip sm V, [by) 



2 4 

2( 

2 V 3 



P,, = -UsmH-l). (70) 



Using (49,53,54) one can obtain: 



dPsiar) ^ -^/^Qj^sin^ T}cos{2n ft + 2ipf ) x [sin 299 cos 2f]t - cos 2(/? sin 217*] (71) 
dt 8 

Prom ()71() it follows that slowly oscillates with the frequency (0,f — 0,). 

In the ideal case, when 0,f = 0, (as it is proposed in jl7j for EDM measurement) (|71() converts to: 

^^^^"^^ = --AnTs'm'^^cos{2nt + 2ipf)sm{2nt-2(p) (72) 
ut 8 

In the general case, when the phases ff and (f are arbitrary, ()72() contains terms that do not depend 
on time and, therefore, P3 linearly grows with time, like the signal from the EDM does. 
It is interesting that making ipf = —ip provides: 

— ^J— ^ ~ cos(2m - 2ip) sin(20t -2ip) = - sin(4m - Aip) (73) 

that makes this contribution to P3 quickly oscillating and depressed. But even in this ideal case it rests 
the contribution caused by the tensor magnetic polarizability ((TJ (in real situation 17 7^ Oj, though). 

Measurement of these contribution provides to measure the tensor electric polarizability. 

According to evaluations |25 ~ 10"^'^ cm^, therefore, for the field -Ee// = B ^ 10^ gauss the 
frequency VIt ~ 10~^ sec~^. When considering modulation one should estimate AQ.T ~ Or(— )^, then 



for the deuteron EDM d = 10~^^ e • cm. 



suppose sec , that exceeds the magnitude of ujd 



6 Particle spin rotation in a storage ring with an internal target 

Let now a gas target be placed inside the storage ring. According to the section 5.1, spin evolution in 
this case is also influenced by the nuclear pseudoelectric and pseudomagnetic fields. 
Time evolution of the spin and tensor polarization can be found by the equations: 

^ (^(t)il>(0) {^{m.u\^{t)) 

Qik is the tensor of rank two (the tensor polarization), for 5 = 1 the tensor = '^{SiSk + SkSi — ^6ik) ■ 
But it should be taken into account that multiple scattering, which occurs in the target, brings to 

both change in the angular spread and depolarization of the particle beam. 

So the spin density matrix formalism should be used to derive equations for evolution of the particle 

spin. 

The density matrix of the system " deuteron+target" is 

P = Pp® Pt, (75) 
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where pd is the density matrix of the particle beam and pt is the density matrix of the target. 
The equation for the particle density matrix can be written as |5l I18j: 



dpp 
dt 



Pp 



+ 



dpp 
dt 



(76) 



col 



where Hp is the Hamiltonian describing particle motion in external electromagnetic fields, 
describes evolution of the density matrix due to collisions with the target atoms (nuclei) . 

can be expressed as follows: 



col 



The collision term ( 



col 



dpd 
dt 



vNSpt 



col 



2m 
T 



[F{9 = 0)p - pF+{9 = 0)] + j dnF{k')p{k')F+{k'] 



(77) 



where k = k -\- q, qis the momentum transferred to a nucleus of matter from the incident particle, v 
is the speed of the incident particles, is the atom density in matter, F is the scattering amplitude 
depending on the spin operators of the deuteron and the matter nucleus (atom), F^ is the Hermitian 
conjugate for the operator F. The first term in ()77() describes coherent elastic scattering of a particle 
by matter nuclei, while the second term is for multiple scattering. Let us consider the first term in 
((77jl. It can be written as: 



dpp 
dt 



2m 



f{o)Pp - PpKoy 



(1) 

= vN- 

col ^ 

The amplitude /(O) of particle scattering in the target at the zero angle is 

/(O) = Sp,F(0)pt. 

As a result one can obtain: 

(1) 

col 



where 



2T:h^ 



eff 



NfiO) 



2TTh^ 



mj rn-7 
The latter expression can be rewritten as: 

27r;i2 



N{d + AiSPt + A2{Sn){nPt) + di{Sn) 



eff 



P' 



-Nd-^SG 



2T:h^ 



NdiiSn) 



(78) 



(79) 



(80) 



(81) 



(82) 



where G 



rwyfi 

to di is absent. 



N{AiPt + A2n{nPt)) (see (j26|) ). For particles with the spin ^ the term proportional 



Finally, the expression (f76|) reads: 



dpp 
dt 



Hp, Pd 



h 



VeffPp - Pp^eff] + ^^^Sp 



j dnF{k')p{k'] 



(83) 



The last term in the above formula, which is proportional to Spf., describes the multiple scattering 
process and spin depolarization aroused from it. 

When studying interaction of particles with a target inside a storage ring the density of the target 
is chosen to make multiple scattering in it small for the observation time; depolarization of the beam 
also appears small. 
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For high-energy particles scattering at smah angles is important. Therefore the large number of 
scattered particles remain on their orbits. For this case the latter term should be considered. 

Further analysis of equations (|8I-{|) in the present paper will be done for such time of the experiment 
that provides the last term in H83() can be cast out. 

As a result we have: 



dpp i 



Hp, Pd 



^{yeffPp-PpK'ff). (84) 



Now the equation for evolution of spin properties of a particle in a storage ring can be obtained. 
The polarization vector is as follows: 

P ^ Sp ppjt) I ^ Sp ppjt) s 

Sppp(t) I{t)S ' ^ ' 

where I{t) = Sp Pp{t) is the beam intensity. 

From l|85p one can get the differential equation providing to find the beam polarization: 

rfP _ Sp ^ 5 1 dljt) 

dt SI{t) I{t) dt ^ ' 

When the particle spin is 5" = 1 then the polarization tensor should also be found: 

p _ Sp pp Qik , . 

Change in the tensor polarization with time can be written as: 



dt I{t) \dt ^'V "'l(t) dt ■ 

Using the equation for the density matrix (|83|) and the expression ()86I88|1 one can obtain the 
equation that summarizes BMT equation (|3U|) for the case when spin evolution is influenced by the 
pseudoelectric and pseudomagnetic nuclear fields. 

Let us consider first particles with the spin 5 = ^. The density matrix for such a particle can be 
expressed as follows: 

Ppi^=n_{t){\i + PS), (89) 

where I is the unit matrix in the spin space. 

With the help of (|85j) it is possible to find that 

= -f>v{yeffPp^^ - Pp^K^ff) = -{^+-^G,P)I.{t), (90) 

where x = vNatot, v is the particle speed, atot is the total cross-section of particle scattering by a 
nonpolarized nucleus and G2 is the imaginary part of the pseudomagnetic nuclear field. The expression 
(POjl can be rewritten as: 

dli{t) 

= ->^abs{P)h{t), (91) 

dt 2 

= X + Pgt, 



where the absorption coefficient x^bs = vN atot + ^(^i{PPt) + ^(^2{P'n){nPt) 
gt = ^vN{aiPt + a2n){nPt). 

The cross-section ai = crj^ti^ Pt) — o'toti''^ -L -Pt) is the difference between the cross-sections of 
particle scattering by a polarized nucleus with the P parallel (P || Pt) and antiparallel (P || Pt) to 
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the target polarization vector Pt in conditions when the particle momentum is orthogonal to the target 



polarization vector (n _L Pt). The cross-section a2 = 0"jJj(n||Pt) — algt{n\\Pt) is the difference between 



the cross-sections of particle scattering by a polarized nucleus with the P parallel (P || Pt) and 
antiparallel {P || Pt) to the target polarization vector Pt in conditions when the particle momentum 
is parallel to the target polarization vector (n||Pt). 

The equation H91|l describes the well known fact that the coefficient of particle beam absorption in 
a polarized target depends on the respective orientations of the beam and target polarization vectors. 

Let us obtain now the equation that describes evolution of the particle vector polarization under 
action of pseudomagnetic nuclear field G. 

According to (|85|86|) one can obtain: 



dP_ 
'dt 



SI 



—Sp^S- 
. (t) dt 



P 



1 



dli{t) 



Ii (t) dt 



2n 



Gi X P 



PiG2P)), 



(92) 



where vecGi = ReG, G2 = ImG. 

The equation ()92() can be rewritten as: 



dP 
It 



PxQr 



{gt - PiPgt)), 



(93) 



where 0,nuc is the frequency of spin precession in the pseudomagnetic nuclear field 



2fx 



Gi 



2'Kh 

rwy 



Ar(Re^iP + ReA2n{nPt)). 



(94) 



Adding the contribution from the pseudomagnetic nuclear field (|93j) to the BMT equation (|30j) one 
can finally obtain the equation describing spin evolution for particles moving in a storage ring with a 
polarized target inside: 



dP 
'dt 



P X {n{d) + ftnu. 

dli{t) 



dt 



- {gt - P{Pgt)), 
-{>c + Pgt)n_{t), 



(95) 
(96) 



where gt = \vN[aiPt + a2n{nPt)). 

The contribution from the pseudomagnetic nuclear field to spin evolution of the 5 = 1 particle 
(for example, deuteron) can be obtained in the similar way. The spin density matrix for the spin 1 
particle is expressed as follows: 



p, = I,{t)i^i+^{PS) + ^P,kQik) 



(97) 



The effective potential energy of deuteron interaction with a polarized target is expressed as 



Nd 



^SG 



NdiiSn)'^ (see ^). 



contains both 



m7 o m-y 

The energy Veff for a deuteron (in contrast to a particle with the spin S = ^ 
deuteron interaction with the pseudomagnetic field and pseudoelectric field (the term proportional to 
di). The equations describing deuteron spin rotation (birefringence effect) in the pseudoelectric field 
of the nonpolarized internal target of a storage ring were obtained in [THl . 

If the internal target in a storage ring is polarized, then using (|82l84l97jl one can obtain the change 
in the beam intensity and polarization vector as follows: 



dhjt) 
dt 



x + 2gtP)-^{2 + Piknink) h{t) 



(98) 
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where x = — ^'^^^ Imdi = —vN{aM=i — o"Af=o), x = vN(Tm=o, o"m=i and (Jm=o are the total cross- 
sections of deuteron scattering by a nonpolarized nucleus for the deuteron state with the magnetic 
quantum number M = 1 and M = 0, respectively (the quantization axis is directed along n). 
The polarization vector can be expressed as: 



dP 

dt 



P X l1„„e] - iPikQtk - Im + 2p{Pgt) + §[n X n']+ ^^^^ 
+f (n(n ■P)+P)- ^P - |(n • n )P 

where rj = -^Re di, n[ = PikUk- 

According to (|99|) the nuclear pseudomagnetic field causes deuteron spin precession with the fre- 
quency ^nuc- 

Contribution H99|) should be added to the equation Ij^Uj) to make them describing evolution of the 
deuteron spin in a storage ring at presence of target. 

The similar approach provides to find the addition to the tensor polarization evolution (see 
expression (|S5|) ). But it is too bulky to adduce it here. 

However, when particle beam absorption in the target can be neglected, one can make all changes 
in the equation, which describes deuteron spin behavior, by replacement of il(ci) with new value 

Note that as the target inside the storage ring is of the finite size then the particle moving in the 
storage ring interacts with the target at times (the expressions for Q-nuc and gt contain the density A^, 
which means the density in the point of particle location i.e. = N{t) ). This is the reason for ^nuc 
and gt oscillating with the frequency of particle rotation in the storage ring Jig. 

For further analysis let us expand VLnuc and gt into Fourier series and be confined with the zero 
harmonics. In this case ^nuc and gt appears to be constants and can be written as: 

^nuc = ^N{ReAiPt + Re^2n(nPt))^ = 

f^f (Re^iPt + Re^2n(nP0) = (100) 



g^(ReAiPi + Re^2n(nPi)), 



m7 L 
m7 

where I is the target length, jt = N ■ I is the target density in cm~^ (usually for a polarized gas target 
jt ~ 10^^ cm~^), T is the period of beam rotation in the storage ring, z/ is the frequency of particle 
rotation in the storage ring, L is the particle orbit length inside the storage ring. 
The zero harmonics for gt is as follows: 

gt = ^vN{aiPt + a2n{nPt))^ = ^^^^^ 
= ijXo-i-P* + cF2n{nPt)). 

Let us evaluate now the effect magnitude. For protons, antiprotons and deuterons with the energy 
from MeV to GeV Re^i ~ lO'^^ ^ iQ-n 

cm. Considering v ~ 10^ s ^ one gets for the nuclear 
precession frequency ^nuc = 10~^ 10~^ s~^. 

Note (see ()16I17() ) that the ratio ^UgAk-^) is proportional to the T-matrix and, as it follows, ^nuc 
depends on the particle energy only due to possible dependence of the T-matrix on energy. 

The obtained estimation for Vlnuc allows to expect, for example for COSY, to observe the spin 
rotation angle "d = ^nuct ^ 10~^ 10~^ rad in the observation time t ~ 10^ s. This value is quite 
observable. 

To prevent suppression of the the spin precession in the pseudomagnetic nuclear field by the storage 
ring magnetic field B the polarized target should be placed in the straight section of the storage ring, 
where the field B = and the particle moves along the straight trajectory. 

Let us now consider the particular example. Suppose the axis OZ is orthogonal to the orbit plane 
{OS\\B) (see Fig.2). The pseudomagnetic field is directed along the axis OX. In this case for the 
particle in the straight section the storage ring the vertical spin component rotates around the direction 
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of the pseudomagnetic field in tlie ZOY plane. Just the change of the vertical spin component should 
be observed. 

As the typical spin rotation angle "d = ^nuct <S 1, then change of the vertical spin component with 
time can be expressed as: 

Psit) = PsiO) - \^''t^ = P3{0) - J^L/- (102) 

The effect can be strengthened when adding a magnetic field b ^ Gi directed along the pseu- 
domagnetic field Gi (along the axis OX in the case under consideration). In this case, the an- 
gle of rotation is expressed as 1) = "dmag + '&nuc = i^magib) + ^nuc)t and the vertical component 
Psit) = PsiO) — \'0'^nag^^ ~ '^mag'&nuc (the terms ~ "^nuc ^^^^ neglected comparing with the previous 
terms). 

When make the direction of the field b changing during the experiment (i.e. the sign of 'dmag 
changing) one gets the possibility to measure the difference ^3(6 jj Gi,t)—Ps{b ][ Gi,t) = 2't}rnag'&nuc- 
This provides to study the possibility of effect observation and measurement of 'dnuc and, therefore, 
to measure the spin-dependent part of the amplitude of elastic coherent zero-angle scattering. 



7 Conclusion 

In the present paper the equations for spin evolution of a particle in a storage ring are obtained 
considering contributions from the tensor electric and magnetic polarizabilities of the particle along 
with the contributions from spin rotation and birefringence effect in polarized matter of a gas target. 

Influence of the tensor electric and magnetic polarizabilities on spin evolution in the resonance 
deuteron EDM experiment is considered in details. 

It is shown that, besides the EDM, the electric and magnetic polarizabilities also contribute to the 
vertical spin component P3. Moreover, the electric polar izability contributes to the P3 component even 
when the deuteron EDM is supposed to be zero and thereby the electric polar izability can imitate the 
EDM contribution. It is shown that unlike the vertical component of the spin P3 the component P33 of 
polarization tensor does not contain contribution from the electric polarizability, whereas contribution 
from the magnetic polarizability reveals only when the deuteron EDM differs from zero. 

It is also shown that when the angle i3 between the spin direction and the vertical axis meets the 

condition sini? = (cos'd = ^J^), the initial value of P33 appears -P33(0) = 0. As a result, the EDM 
contribution to the measured signal linearly growth in time starting from zero that is important for 
measurements. 

Therefore, measurement of the P33 component of the deuteron tensor polarization seems to be of 
particular interest, especially because appearance of the nonzero component P33 on its own indicates 
the EDM presence (in contrast to the P3 component, which appearance can be aroused by the tensor 
electric polarizability, rather than EDM). 

It is also shown that study of spin rotation and the birefringence effect for a particle in a high 
energy storage ring provides for measurement of the spin-dependent real part of the coherent elastic 
zero-angle scattering amplitude. 



References 

[1] C.Lechanoine-Lelue and F.Lehar Rev.Mod.Phys. 65 (1993) 47. 

[2] R.Engels, K.Grigoriev, A.Gussen et al. Proceedings of the 6th International Conference on 
Nuclear Physics at Storage Rings STORI'OS, Schriften des Forschungszentrums Julich, Mat- 
ter and Materials, Volume 30, 2005, p. 381. 



19 



[3] V.G. Baryshevsky, Sov. J. Nucl.Phys. 38 (1983) 569; V.G. Baryshevsky, Phys. Lett. B 120 
(1983) 267. 

[4] V.G. Baryshevsky, I.Ya.Dubovskaya, Phys. Lett. B 256 (1989) 529. 

[5] V.G. Baryshevsky, A.G. Shekhtman Phys. Rev. C 53, n.l (1996) 267. 

[6] V. G. Baryshevsky, Phys. Lett. 171A (1992) 431. 

[7] V. G. Baryshevsky, J. Phys.G 19 (1993) 273. 

[8] V. G. Baryshevsky, K. G. Batrakov and S. Cherkas J. Phys.G 24 (1998) 2049. 

[9] V. Baryshevsky, Proceedings of the 6th International Conference on Nuclear Physics at 
Storage Rings STORr05, Schriften des Forschungszentrums Julich, Matter and Materials, 
Volume 30, 2005, p.277. 

[10] V.G. Baryshevsky and M.I. Podgoretsky, Zh. Eksp. Teor. Fiz. 47 (1964) 1050. 

[11] A. Abragam et al., C.R. Acad. Sci. 274 (1972) 423. 

[12] M. Forte, Nuovo Cimento A 18 (1973) 727. 

[13] A. Abragam and M. Goldman, Nuclear magnetism: order and disorder (Oxford Univ. Press, 
Oxford, 1982). 

[14] W. Bryan, Montague Phys. Reports 113, no.l (1984) 1-96. 

[15] S.R. Mane, Yu.M. Shatunov and K. Yokoya, Rep.Prog.Phys. 68 (2005) p. 1997-2265. 
[16] F. Farley et al., Phys. Rev. Lett., 93 (2004) N 5. 

[17] Yu. Orlov, Proceedings of the 6th International Conference on Nuclear Physics at Storage 
Rings STORr05, Schriften des Forschungszentrums Julich, Matter and Materials, Volume 
30, 2005, p. 223. 

[18] V. Baryshevsky, A. Shirvel, LANL e-print archive: ||hep-ph/0503214 

[19] V. Baryshevsky, LANL e-print archive: (hep-ph / 0504064[ 

[20] V. Baryshevsky, A. Gurinovich LANL e-print archive, h ep-ph/ 0506135] 

[21] Jiunn-Wei Chen, Harald W. GrieBhammer, Martin J. Savage and Roxanne P. Springer, 
Nucl.Phys. A644 (1998) 221-234 (LANL e-print arXiv:nucl-th/9806080J . 

[22] M.Lax, Rev.Mod.Phys. 23 (1951) 287. 

[23] M.L. Goldberger and K.M. Watson, Colhsion Theory: John Wiley, New York, 1964. 



20 



